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ABSTRACT 

We have used the He I discontinuities at 3421 A to determine the electron 
temperatures, designated T e (He i), for a sample of five Galactic planetary nebu- 
lae (PNe). We compared T e (He i) with the electron temperatures derived from 
the hydrogen Balmer jump at 3646 A, designated T e (H i), and found that T e (He i) 
are generally lower than T e (H i). There are two possible interpretations, a) the 
presence of substantial He 2+ zone, or b) the presence of hydrogen-deficient cold 
clumps within diffuse nebulae. A series of photoionization models were con- 
structed to test the two scenarios. We found that the observed T e (He i)/T e (H i) 
discrepancies are beyond the predictions of chemically homogeneous models. Our 
modelling shows that the presence of a small amount of hydrogen- deficient in- 
clusions seems to be able to reproduce the observed intensities of He I discon- 
tinuities. We stress the value of He I discontinuities in investigating nebular 
physical conditions. Albeit with some observational and technical limitations, 
He I discontinuities should be considered in future modelling work. 

Subject headings: ISM: general — planetary nebulae: general 

1. Introduction 

The accurate determination of chemical abundances of planetary nebulae (PNe) is of fun- 
damental importance to the understanding of the nucleosynthesis of low- and intermediate- 
mass stars and the chemical evolution of galaxies. However, one of the main problems in 
nebular astrophysics is that the heavy element abundances derived from collisionally excited 
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lines (CEL) are often lower than those derived from optical recombination lines (ORLs). 
The typical ORL/CEL abundance discrepancy factor (ADF) is ~ 2 for general P Ne. The 
most extreme case to date is Hf 2-2, which has a ADF of abou t 70 (ILiu et al.ll20061). Recent 
reviews on this problem have been presented by lLiul (120061 ) and lPeimbert &: Peimbertl (120061 ) 
(see also the references therein). There are two possible solutions: (1) the presence of tem- 
perature and density variations in chemically homogeneous nebulae; (2) the two-component 
nebular model with hydrogen- deficient inclusions embedded in the diffuse nebula. In order 
to test the two scenarios, we need to further investigate nebular physical conditions, in par- 
ticular for the regions that ORLs originate in. For this purpose, new plasma diagnostic tools 
rather than the classical CEL diagnostics are badly required. 



Zhang et al.l (120041 ) applied the hydrogen Balmer jump to determine the electron tem- 
peratures - hereafter T e (H i) - for a large sample of PNe. They found that T e (H i) is 
systematically lower than that derived from the ratio of nebular to auroral line s of [O ill] 
- hereafter T e (0 ill), consistent with previous results by ILiu fe Danzigerl (11993). Th e dis- 
crepancies between T e (H i) and T e (0 ill) were first studied by iPeimbertl (119671 . Il97ll ). who 
found that temperature variations within nebulae may lead to higher T e (0 ill) compared 
to T (H i). To quantitatively study the problem, they defined the mean square tempera- 
ture variation, t 2 . Photoionization models of chemic ally and spatially homogeneo us nebulae 
yielded typical valu es of t 2 between 0.003- 0.01 5 (see IPeimbert fe Peimbertl 120061 ). However, 
the observations by lLiu fc Danzigerl (119931 ) and lZhang et al.l (120041 ) indicated to considerably 
large temperature variations, which are far beyond the predictions of typical photoionization 
models. This can be ascribed to either additional energy inputs for chemically homogeneous 
nebulae or chemical inhomogeneities. 



Zhang et al.l (l2005aU bl) developed a method to use the He I recombination line ratios to 
diagnose electron temperatures - hereafter T e (He i). The comparison between T e (He i) and 
T C (H i) is essential to discriminate between the two- abundance model and chemically homo- 
geneous nebulae with temperature and density variations as responsible for the CEL/ORL 
abundance problem since the two scenarios p redict different relat ions between T e (He i) and 
T e (H i). From a study of a sample of 48 PNe, Zhang et al. (j2005a ) (hereafter Paper i) found 
that T e (He i) is significantly lower than T P (H i), in favor of the t wo-ab undance model. How- 
ever, there are two arguable problems. IPeimbert Sz Peimbert found that for some 
PNe the T e (He i) values derived from the He I AA3889, 4471, and 7069 are different with 
those derived from the He I AA6678,7281. This might be partially caused by the uncer- 
tainties of the density of H e + zones and/or the optical depth effects on the He I A3889. 
Peimbert fc Peimbertl (120061 ) proposed that at least 10 different He I lines are required to 
accurately deduce these parameters. Moreover, for high-excitation PNe the He + and H + 
zones are not identical, and the presence of a substantial He 2+ zone may partially contribute 
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to the difference between T e (He i) and T e (H i) 



( Garcia- Roj as et al. 


2005; 


Guseva et al. 


2006, 


2007) 



temperature differences in H II regions are not so obvious as those found in PNe. 

In this paper, we use an alternative method, fitting the He I discontinuities, to determine 
T c (He i) for five Galactic PNe. This method is made possible by the availability of high signal- 
to- noise spectra of these nebulae resulting in a clear detection of the helium jump. The 
current study allows us to further investigate the relation between T e (He i) and T e (H i). We 
examine the two-abundance models and photoionization models of chemically homogeneous 
nebulae as the explanation of the observed T e (H i)/ T c (He i) discrepancies. This paper is a 
comple ment to the authors' work on recombination spec tra as probes of physical conditions 
of PNe f lLiu fc Danziger1ll993l : Izhang et al.lbooi l2005aB ). 



2. Observations and data reduction 

The sample includes two low-excitation PNe and three high-excitation PNe (see Tabled]). 
The long-slit spectra were obtained using the f-18 Nasmyth focus of the 2.3 m Advanced 
Technology Telescope at Siding Spring Observatory on March 2000. Both arms of the Double 
Beam Spectrograph (DBS) were used with the 1200B grating in the second order to give 
spectral resolution of 0.3 A/pixel. The dichroic may be removed or replaced with a flat mirror 
allowing the use of either arm independently. The detector was a SITE 1752 x 532 chip, 
giving a spectral coverage from 3240-8520 A. In order to avoid loss of light due to atmospheric 
refraction each PN was observed with the spectrograph slit aligned to the parallactic angle. 
Furthermore, all the observations were carried out under the best seeing conditions, along 
with airmass around 1.2 to ensure the total UV light into the spectrograph slit. For each 
PN, the exposure time was 1500 s. 

The LONG92 package in MIDAS0 was used for the data reduction following the standard 
steps, including bias subtraction, flat-field correction, cosmic-ray rejection, and wavelength 
calibration with a helium-argon calibration lamp. Flux calibration was conducted using 
the spectroscopic standards EG 131 and LTT4364. Finally, the spectra were dereddened 
using the stand ard Galactic extinction law for a total-to -selective ext i nctio n ratio of R = 3.1 



( ]Howarthlll983l ) and the extinction coefficients given by ICahn et al.l (Il992l ). 



midas is developed and distributed by the European Southern Observatory. 
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3. Analysis and results 



The spectra allow us to simultaneously determinate T e (He i) and T C (H i) in terms of 
the He I discontinuities at 3421 A (Jn e = /3421- — ^3421+), produced by He + recombination 
to the He I 2p 3 P° level, and the H I Balmer discontinuities at 3646 A (J H = ^3646- ~~ ^3646+)- 
Note that we do not distinguish between the electron temperatures derived from the He I 
discontinuities and from the He I lines and u se TpjHe 1) to denot e both of them. The method 



is similar to that used to obtain T e (H 1) in IZhang et al.l (120041 ) . To synthesize the spectra 



we consider the contributions from the free-free, free-bound, and bound-bound emission 
from H I, He I, and He II and the two-photon decay from the 2 2 S level of hydrogen. The 
coefficients of the H I, He I, an d He II continuous emission as functions of temperature are 
available in lErcolano fc Storey! (120061 ) . The scattering from the central star is simulated 
using a modified black-body function, I\ oc A~^I?(A, T), where (3 is the spectral index {(3 = 
for stardand black-body spectrum), B(\, T) is the Planck fun ction, and T is the temperature 
of the central star. We refer the reader to lZhang et al.l (120041 ) for a detailed description of our 
fitting procedure. In order to illustrate our method, Fig. [T] gives the theoretical spectra at 
different T e (He 1). It is clear that increasing of T e (He 1) will decrease the He I discontinuities. 
This is a consequence of the fact that the number of the lowest-energy free electrons, which 
is proportional to the intensity of the He I free-bound emission at 342 1~ A, decreases with 
increasing electron temperature. For the same reason, the H I Balmer discontinuities decrease 
with increasing T e (H 1). Note that the Jh values may slightly change with T e (He 1) because 
of a small contribution from the He I discontinuity at 3680 A produced by recombination to 
the 2 X P level. 

Consequently, T c (He 1) and T C (H 1) can be derived by fitting the theoretical to the 
observed spectra. Before the fitting, we require to measure the He + /H + and He 2+ /He + 
abundance ratios. The He + /H + ratios are obtained from a few strong He I lines around the 
He I discontinuities. The deduced He + /H + ratios are relatively insensitive to the adopted 
electron temperatures and densities. The fl ux uncertain ti es of the He I lines are less than 
15%. The He 2+ /He + ratios are taken from Cahn et al. ( 19921 ) (see Table [T]). The errors 
in temperature determination caused by uncertainties of the He 2+ /He + ratios are negligible 
since the contribution from He II to the continuum emission is much smaller than that from 
H I and He I. 

Figure. [2] shows the fitting results for the five PNe. In order to reduce the errors caused 
by uncertainties in the reddening correction and flux calibration, we have normalized all the 
spectra to H 11 A3770. The resultant T e (He 1) and T e (H 1) are given in Table [TJ Given 
the relations T e (He I, H I) oc Jh^h", where a ~ 1.5, the errors introduced by measurement 



- 5 - 



uncertainties can be estimated through 

AT e (He I, H I) = — ^— T c (He I, H I), (1) 

^He,H 

where a is the rms noise derived from the spectral regions free of emission lines, which cover 
a wide wavelength range. We found that for all the nebulae the 3a values are well below the 
</He,H values. Another error source arises in the uncertainties of the He I line emission coeffi- 
cients, which might introduce an error in the He I temperature determin ation of about 20% 



according to a comparison of the atomic data gi ven by different authors (IBrocklehurst et al. 



19991 ; iBenjamin et al.lll999l ; iBauman et al.ll2005l ). For high-excitation PNe, additional errors 



are caused by the strong O ill Bowen lines around the He I discontinuities. In such cases, our 
fittings were based on the spectral regions with wavelengths departing from 3421 A, intro- 
ducing a ~ 15% errors in derived T e (He i) caused by uncertainties in flux calibration. The 
emission efficient of hydrogen two-photon continuity was deduced under Case A assumption 
and might be significantly underestimated under Case B conditions due to the collisional 
excitation from 2 2 P to 2 2 S. However, the excess emission of two-photon process can be 
mostly compensated by the simulated scattering continuities, and we do not require accu- 
rate estimate of the absolution magnitude of two-photo emission. The uncertainties in the 
smooth continuities hardly affect our results mainly based on the 3421 A discontinuity. 

The right panels of Figure. [2] give zoom-in views of the spectra around the He I discon- 
tinuities. For some PNe, the narrow spectral ranges from 3421-3500 A cannot be well fitted 
by our models due to the facts that the models did not treat the extended wings of O ill 
Bowen lines and the blended high- J He I lines at the red side of the 3421 A discontinuity, and 
the low sensitivity near the CCD edges leads to a larger slope. Nevertheless, these effects on 
our results should be somewhat insignificant since the deviations exist only in a relatively 
limited wavelength range. Of course, if our models merely take into account the spectra 
from 3300-3600 A, we can achieve a better fitting in such a narrow wavelength range by 
adjusting the contribution from scattering continuities. However, this will result in a failed 
fitting of the global spectra. Modelling different spectral ranges to a reasonably good extent, 
we have assessed the errors in T e (He i). In Figure. [21 we overlaid the best fitting to the nar- 
row spectral ranges (3300-3600 A) of NGC2440, NGC3242, and NGC3918, which results in 
T e (He l)= 6840 K, 7500 K, and 8250 K for the three high-excitation PNe, respectively. We 
noted that these values, which should be larger than the real T e (He i) values, are still lower 
than T e (H i). 

Three of the PNe in the recent sample have been studied in Paper I. For comparison, 
Table [1] gives T e (H i) and T e (He i) determined in Paper I. Considering the error bars, they 
are in agreement with the results derived in the recent work. Note that the T e (He i) in 
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Paper I were derived from He I lines ratios and might possibly be different with those from 
He I discontinuities due to the effects of radiative transfer and temperature inhomogeneities. 

In Fig. [31 we compare T e (H i) and T c (He i). These points derived in the recent work and 
those derived in Paper I are located at approximately the same position in this figure. An 
inspection of Fig. [3] shows that T e (He i) is consistently lower than T e (H i), confirming our 
previous results in Paper I. In the scenario of the two-abundance model, the He I emission are 
significantly contaminated by emission from the extremely cold hydrogen-deficient inclusions 
while the H I recombination spectra dominantly originate in the normal component under 
a typical temperature of T e ~ 10 4 K. This provides a plausible explanation for the lower 
r fi (He i) compared to TJH i) . In Fig. [31 we also overplot the results for H II regions obtained 



by iGarcfa-Rojas et al.l (120051 ). which are located at a more compact area in this figure and 



indicate to a smaller T e (He l)/T c (H i) discrepancy in H II regions compared to PNe. 



4. Discussion 
4.1. He I lines 

A few strong He I lines have been clearly detected in the red arm. They allow an 
independent examination of our results derived from the He I discontinuities. This is non 
trivial in that density inhomogeneities and radiative transfer effects could affect our results 
to some extent. 

Here we took into account the four strongest He I lines in the red arm, AA7281, 6678, 
7065, and 5876, as shown in Fig. HI Their intensities are listed in Table [2j In Paper I, we 
have shown t hat the He I ratio I(A72 81)/I(A6678) provides an useful temperature diagnostic. 



According to lBenjamin et al.l (120021 ). the He I AA7065, 5876 lines are from triplet states and 
are essentially affected by the optical depth of the 2s 3 S level. Fig. gives the theoretical 
predictions of the He I ratio J(A7281)//(A6678) versus He I ratio J(A7065)//(A5876) at 
different electron temperatures, electron densities (N e ) , and optical depths of the He I 2s 3 S- 
3p 3 P°A3889 line (T3889). The observed values for our sample are also shown in this figure. 

The diagram is useful to investigate nebular physical conditions. An inspection of 
Fig. [5] shows that the T e (He 1) values are consistently lower than 10 000 K, which represents 
the lowest value of T e (H 1) in our sample, supporting the conclusion drawn from the He I 
discontinuities. We also found that the r 3889 values range from to 10. These line ratios have 
only weak dependence on electron densities and are hardly used to determined iV e (He 1). 

Table [2] gives the T e (He 1) values derived from the He I line observations. The uncer- 
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tainties of electron densities and optical opacities introduce an error of about 1500 K. Using 
the spectral lines falling in the overlapping wavelength region of two adjacent orders, we 
estimated the uncertainties in flux calibration to be about 5%. The combined errors of the 
resultant T e (He i) amounts to about 2200 K. The T e (He i) deduced in this paper are biased 
to larger values compared to those in Paper I (although they are marginally consistent with 
each other considering the error bars). This is partly attributed to the fact that in the 
present observation the slit was put in a position closer to the center. However, the primary 
cause might be that the uncertainties in flux calibration of the red arm are probably larger 
than expect. We performed a comparison of the present spectra and our previous long-slit 
spectra and found that the former have systematical higher line fluxes in very red wave- 
length regions. Therefore, we suppose that the fluxes of the He I A7281 might have been 
significantly overestimated, and the T e (He i) in Table [2] should be treated as upper limits, 
which are still lower than T e (H i). A detailed discussion of He I lines requires more precise 
observations in the red arm and is beyond the focus of the present paper. 

One should bear in mind that the electron temperatures derived from the He I disconti- 
nuities and those from the He I line ratios are not necessary identical since they have different 
temperature dependences, and thus weight different regions if large temperature fluctuations 
exist. The comparison of He I discontinuity- and line-temperatures is a potentially intriguing 
subject for future research. 



4.2. Model tests 

An issue is invited, whether the large T e (H i)/ T e (He i) discrepancies in high-excitation 
PNe can be ascribed to the presence of substantial He 2+ zone which has a higher temper- 
ature than He + zone. In order to investigate the problem, we first compare the observed 
results with the predictions by ID photoionization models of chemically homogeneous neb- 
ulae. We then construct 3D photoionization models to test the two-abundance model as the 
explanation of the temperature discrepancies. 



Jf.,2.1. Chemically homogeneous model 



To date, the published photoionization models for certain PNe with a chemically homo- 
geneous medium predict only minor discrepancy between T e ( He i) and T e (H i). For ex ample, 
the classic photoionization model of NGC 7662 constructed by Harrington et al.l (119821 ) yields 
T C (H i)— T e (He l)~ 600 K. To investigate the case of general PNe, we use the cloudy code 
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( IFerland et al.l Il998l ) to construct a series of ID photoionization modellings of chemically 
homogeneous nebulae. In order to compare the modelling results with the observed ones, we 
define modelled T e (He i) and T e (H i) by 



T e (He I, HI) 



/ T e (r) Jue,u{r)dr 

(2) 



j T c (r)N c (r)N Hc+!il+ (r)T c (r)~ 2 / 3 dr 



JiV c (r)iV Hc+iH+ (r)r c (r)-2/3rf r ' 

where r is the position vector, N e (r) is the electron density, and iV He + H +(r) is the number 
density of He + or H + atoms. Then T e (He i) and T e (H i) can be simulated from the modelled 
T e (r), N e (r), and A^ He + H +(r). The temperatures defined in Eq. (2) are similar to those 
derived from the observed He I and H I discontinuities. 

For the modellings, a black-body spectral energy distribution with a luminosity = 
10 38 ergs -1 is assumed for the central star. We assume that the nebula has a homogeneous 
density structure with a hydrogen density of n(H) = 10 3 cm -3 . The assumed n(H) represents 
the lower limit of typical PNe. At a given UV radiation field, increasing n(H) will decrease 
the ionization parameter and the He 2+ /He + ratio, and thus decrease the contribution of 
He 2+ zone to the T e (He l)/T e (H i) discrepancies. 

Our modellings reproduce the T e (He i) and T e (H i) values in different nebular metallic- 
ities (Z) and effective temperature of the ionizing star (T*). The results are illustrated in 
Fig. [31 These modellings show that the presence of He 2+ zone can cause a higher T e (H i) 
compared to T e (He i). The T e (H l)/T e (He i) discrepancies sharply increase with increasing 
T* and have only weak dependency on Z. Nevertheless, the large differences between the 
observed T e (He i) and T C (H i) are obviously beyond the model predictions although the claim 
is probably not applied for H II regions. An inspection of Fig. [3] suggests that chemically 
homogeneous models are unable to reproduce the extremely low T e (He i) measured in PNe. 
The observed T e (H i) values can be well reproduced by using appropriate T* and Z. How- 
ever, in order to match the observed T e (He i) for most PNe, we require extremely high T* 
(S> 200, 000 K) and Z (^> 1.5 Z & ), which are unlikely real for our sample. 

Our models do not include the photoelectric heating by dust grains , which can efficiently 



enha nce the temperature in the nebular regions close to the central star (jStasihska fc Szczerba 



200ll ). In addition, there might be other unexpected heating sources that lead to greater tem- 
perature fluctuations. It is too early to conclude with certainty that chemical inhomogeneities 
must be introduced to explain the observed temperature discrepancies. In any case, a more 
comprehensive photoionization model including extra heating mechanisms should take into 
account the He I discontinuities. 
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4-2.2. The two-abundance model 
To test the two-a bundance model, we use the 3D photoionization code mocassin 



( lErcolano et al.ll2003al ). which is capable of treating nebular geometric asymmetries, den- 
sity and composition inhomogeneities as well as the diffuse radiation fields self-consistently, 
and thus enables to more rigorously study the physical properties and chemical composition 
of the postulated hydrogen-deficient component compared to ID model. A series of two- 
abundance models are constructed and the reproduced Jr and Jn e are compared with the 
observed values. We found that the presence of a very small amount of hydrogen-deficient 
material is enough to interpret the observations. An example is given in Fig. 0, where we 
compare the normalized Jh and Jh e observed for our PN sample and those predicted by 
one of the mocassin models. For this two-abundance model, we have assumed that the 
diffuse nebula has a spherical structure, a n.(H) value of 1 3 cm" 3 , and a metallicity of Zpn, 



where Zpn is the average metallicity of the Galactic PNe (IKingsburgh fc Barlowlll994j ). The 
hydrogen-deficient inclusions have a filling factor of 0.18% and a n(H) value of 70 cm -3 , and 
are enhanced by a factor of 200 in He and by a factor of 800 in CNONe compared to the dif- 
fuse gas. An effective temperature of 120,000 K is assumed for the central star. According to 
our modelling, the electron temperature of the hydrogen- deficient inclusions is about 500 K, 
which is much lower than the average temperature of the diffuse gas of 10,770 K. The cold 
hydrogen- deficient inclusions contribute about 45% and 0.2% of He I and H I line emission 
intensities, respectively. As shown in Fig. EJ the model result can closely match the observed 
J H and J He - 

For comparison, we construct a series of chemically homogeneous models using mo- 
cassin in the parameter ranges 50,000K< T* < 200, 000 K, 500cm~ 3 < n(H) < 8000 cm" 3 , 
and 0.6Z PN < Z < 1.6Z PN . The predicted region of Jh and Jn e is given in Fig. Again it 
is confirmed that even for high-excitation PNe chemically homogeneous models are unable 
to reproduce the observed He I discontinuities. 

The two- abundance model can provide a natural solution for the observed T e (He i) /T e (H i) 
discrepancies. Nevertheless, we have noted some problems in our two-abundance modelling. 
In order to explain the low T e (He i), we require to assume a low density and high enhance- 
ment factor of heavy elements for the hydrogen-deficient inclusions to get effective cooling. 
Such a low-density and low-temperature component is unstable inside the nebulae. More- 
over, our models yield an ADF of 10-20, higher than the typical value for general PNe. 
Probably, these imply to the presence of some additional coolants inside hydrogen-deficient 
inclusions. 



The hydrogen-deficient inclusions postulated in the current study have similar proper- 
ties to Abell 30's knots, which show extremely He-rich (IWesson et al.ll2003l ; lErcolano et al. 
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2003bl ). However, the properties of the hydrogen-deficient inclusions in different PNe are 



probably not identical. A complete understanding of the hydrogen-deficient component re- 
quires detailed modelling and precise measurements of heavy-element ORLs for individual 
PNe, which is beyond the scope of this paper. Our results show that the observations of 
He I discontinuities can provide some constraints on future study of the hydrogen-deficient 
inclusions. 



5. Conclusions 

Based on a high signal-to-noise spectroscopy of a sample of PNe, we explored the pos- 
sibility of using the He I discontinuities at 3646 A as probe of nebular physical conditions. 
T c (He i) is found to be systematically lower than T e (H i), confirming the results of Paper I, 
in which we determined T e (He i) using the He I recombination lines. We have investigated 
the contribution of He 2+ zone on the T e (He i)/T e (H i) discrepancies by constructing a se- 
ries of photoionization models. We have shown that the chemically homogeneous models 
predict too small temperature discrepancies as compared to the observations. The observa- 
tions, however, can be naturally interpreted if one assumes the presence of a small amount 
of hydrogen-deficient cold inclusions within the diffuse nebula. Our results thus provide a 
possible evidence in favor of the two-abundance model as the explanation of the CEL / ORL 
abundance problem. Nevertheless, we cannot completely rule out the scenario of temperature 
fluctuations if some extra heating mechanisms exist. 

To completely understand the properties of the postulated small hydrogen-deficient 
inclusions, we need to construct detailed 3D photoionization models with very fine mesh 
grids and take account of the ORLs from helium and heavy elements. For that, an enormous 
amount of computer time is required. The detailed two-abundance models for individual 
objects will be reported in forthcoming papers (e.g. Yuan et al. in preparation). 

In the current paper, we show that the He I discontinuities, notwithstanding some ob- 
servational and technical limitations, provide a valuable probe of nebular physical conditions 
and can place a tight constraint for the modelling. The method proposed in this paper is par- 
ticularity appropriate for studying low-excitation PNe and H II regions that are free from the 
effects of strong O ill Bowen lines. Instrument development and atomic data improvement 
will be invaluable to provide more robust conclusion. 

The authors are grateful to M. A. Dopita and B. Rocca-Volmerange for their help 
in ensuring successful observations. We thank the anonymous referee for useful comments. 
CTH is grateful to HKU for having made possible the present output. The observations were 
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Fig. 1. — Theoretical H I and He I recombination spectra at different T e (He i). For the 
calculations, we have assumed T C (H l)=10 000K, He+/H+=0.f, and He 2+ /He+=0. 
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Fig. 2. — Left panels: a comparison of observed (solid lines) and theoretical (dashed lines) 
spectra. The dotted-dashed lines and the dotted lines represent the theoretical nebular 
continuum and scattered stellar light, respectively. Right panels: zoom-in views of the 
spectra around the He I discontinuities at 3421 A. The contributions from scattered-light 
continuum have been subtracted. The modelling nebular continuities are overlaid with thick 
solid lines. Note that although our models give the best fitting to the global spectra from 
3240-4000 A, they may not satisfactorily reproduce the narrow spectral ranges of the three 
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Fig. 3. — T e (H i) versus T e (He i). The filled circles with error bars are those with T e (He i) 
determined from the He I discontinuities in this work. The open circles are those with 
T c (He i) deduced fr om He I lines in Paper I. The filled triangles are the results for H II 
regions obtained by iGarcia-Rojas et al.l (120051 ). These lines are the predictions by cloudy 
for different effective temperatures of the central stars and metallicities of the nebulae (see 
the context). 
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Fig. 4.— He I 4471, 5876, 6678, 7281 lines (from left to right) detected in IC 2501, IC 4634, 
NGC 2240, NGC 3242, and NGC 3918 (from top to bottom). 
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Fig. 5.— Theoretical predictions of the He I line ratios 7(A7281)//(A6678) versus 
J(A7065)//(A5876), derived at different physical conditions: T c = 5000 K (squares), fOOOOK 
(triangles), and 20 000, K (stars); N c = 10 4 cm" 3 (light lines) and 1 6 cm -3 (thick lines); 
^3889 = (solid lines) and fO (dashed lines). The filled circles represent the observations. 
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Fig. 6. — Jh//(H11) versus Jn c //(He I A3634). The filled circles denote the observed values. 
The open triangle is the prediction of a two- abundance model (see the text). The enclosed 
region represents the zone occupied by the predictions of chemically homogeneous models 
in the parameter ranges 50,000K< T* < 200, 000 K, 500cm~ 3 < n(H) < 8000 cm" 3 , and 

O.GZpjy < Z < l.GZpjy. 
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Table 1: T e (He i) and T e (H i) 



Source 


He 2+ /He +a 


T c (He l) b (K) 


T.(Hi 


) (K) 






This work 


Paper I 


This work 


Paper I 


IC 2501 


0.002 


4700±3200 




10000±600 




IC 4634 


0.000 


3800±2800 


5000±1000 


10000±500 


8500±400 


NGC 2440 


0.401 


5700±4100 c 




15000±800 




NGC 3242 


0.230 


5000±2800 d 


4500±1000 


12000±1200 


10200±1000 


NGC 3918 


0.303 


5500±3900 e 


6000±1000 


11000±300 


12300±1000 



a From lCahn et al.l (Il992h . 

b Note that the T e (He i) values given by Paper I are derived from He I line ratios. 

c If merely the spectral range from 3300-3600 A is taken into account, we derived a value of 
6840 K, which represents the upper limit of T e (He i) (see the dotted lines in the rigth 
panel of Figure [2]). 

The upper limit is 7500 K. Details indicated in Footnote c. 

e The upper limit is 8250 K. Details are as indicated in Footnote c. 



Table 2: He I line intensities relative to /(H/5) = 1 and T e (He i) derived from He I lines. 



Transition 


IC 2501 


IC 4634 


NGC 2440 


NGC 3242 


NGC 3918 


Singlets 












A6678 


0.0400 


0.0590 


0.0244 


0.0381 


0.0331 


A7281 


0.0086 


0.0142 


0.0056 


0.0082 


0.0068 


Triplets 












A5876 


0.1790 


0.2125 


0.1039 


0.1210 


0.1367 


A7065 


0.0904 


0.1057 


0.0393 


0.0494 


0.0513 


T c (He l) a 


7500 K 


8400 K 


8000 K 


7400 K 


7100 K 



a Bears an error of ~ 2200 K. 



